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Core-collapse supernovae are among Nature's grandest explosions. They are powered by 
the energy released in gravitational collapse and include a rich set of physical phenom- 
ena involving all fundamental forces and many branches of physics and astrophysics. We 
summarize the current state of core-collapse supernova theory and discuss the current set 
of candidate explosion mechanisms under scrutiny as core-collapse supernova modeling is 
moving towards self-consistent three-dimensional simulations. Recent work in nuclear the- 
ory and neutron star mass and radius measurements are providing new constraints for the 
nuclear equation of state. We discuss these new developments and their impact on core- 
collapse supernova modeling. Neutrino-neutrino forward scattering in the central regions 
of core-collapse supernovae can lead to collective neutrino flavor oscillations that result in 
swaps of electron and heavy-lepton neutrino spectra. We review the rapid progress that 
is being made in understanding these collective oscillations and their potential impact on 
the core-collapse supernova explosion mechanism. 

1 Overview: Core-Collapse Supernova Theory 

The ultimate goal of core-collapse supernova theory is to understand the mechanism driving 
supernova explosions in massive stars, connect initial conditions to the final outcome of collapse, 
and make falsifiable predictions of observable signals and explosion features. These include neu- 
trino, gravitational wave, and electromagnetic signals, nucleosynthetic yields, compact remnant 
masses, explosion morphologies, and pulsar kicks, spins, and magnetic fields. 

Baade and Zwicky, in their seminal 1934 article [1 , first hypothesized that a "supernova 
represents the transition of an ordinary star into a neutron star, consisting mainly of neutrons." 
This basic picture still holds today and the road to its refinement has been, at best, meandering 
and bumpy. When the nuclear fuel at the core of a massive star is exhausted, the core becomes 
electron degenerate and, upon reaching its effective Chandrasekhar mass, undergoes dynamical 
collapse. Electron capture on free protons and protons bound in heavy nuclei reduces the 
electron fraction (Yf,; the number of electrons per baryon) and accelerates the collapse of the 
inner core. When the latter reaches nuclear density, p^uc ~ 2.7 x 10^'^gcm^'^, the nuclear 
equation of state (EOS) stiffen^ leading to core bounce and the formation of the bounce shock 
at the interface of inner and outer core. The shock initially rapidly propagates out in radius 

^The stiffening of the EOS near pnuc is due to the repulsive effect of the strong force at small distances and 
r = dlnP/dlnp jumps from ~4/3 to > 2. Neutron degeneracy, which is non-relativistic at bounce, only gives 
r !ii 5/3. 
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and mass coordinate, but the work done to break up infalling heavy nuclei and energy losses 
to neutrinos quickly sap its might. The shock stalls within tens of milliseconds of bounce and 
turns into an accretion shock at a radius of ^100 — 200 km [5 . 

Core collapse liberates ~3 x 10^^ erg = 300Bei/ie of gravitational binding energy of the 
neutron star, ~99% of which is radiated in neutrinos over tens of seconds. The supernova 
mechanism must revive the stalled shock and convert ^1% of the available energy into energy 
of the explosion, which must happen within less than ~0.5 — Is of core bounce in order to 
produce a typical core-collapse supernova explosion and leave behind a neutron star with the 
canonical neutron star gravitational mass of ^IAMq [HH]. 

The neutrino mechanism [21[5] for core-collapse supernova explosions relies on the deposition 
of net neutrino energy (heating > cooling) in the region immediately behind the stalled shock, 
heating this region and eventually leading to explosion (for details, see the excellent discussion 
in [H]). While having great appeal and being most straightforward, given the huge release of 
neutrino energy in core collapse, the simplest, spherically-symmetric form of this mechanism 
fails to revive the shock in all but the lowest-mass massive stars (0-Ne cores) [71-[TUj. 

Indications are strong that multi-dimensional effects, principally turbulent convective over- 
turn and the standing-accretion-shock instability (SASI, e.g., [HKT^] and references therein) in- 
crease the efficacy of the neutrino mechanism by boosting neutrino heating [I3HI6] or, as suggest 
by [5], by reducing neutrino cooling. This is generally borne out by recent fully self-consistent 
axisymmetric (2D) neutrino radiation-hydrodynamics simulations with an energy-dependent 
treatment of neutrinos, but their detailed results vary significantly from group to group and a 
clear picture has yet to emerge. Marek et al. [T7j reported the onset of explosions in a nonrotat- 
ing 11.2-M0 (at zero-age main sequence [ZAMS]) and in a slowly spinning IS-M© star, setting 
in at ~200 ms and ~600 ms after bounce, respectively, and the estimate explosion energies 
are on the lower side of what is expected from observations. The exploding simulations used 
the softest variant of the EOS by Lattimer & Swesty (LS) [TS] and included a stronger quasi- 
relativistic monopole term in the gravitational potential. A similar, so far unpublished |19j . 
calculation of core collapse in a 11.2-M0 star with the stiffcr EOS of H. Shen et al. |20] also 
produced an explosion while simulations with the very stiff EOS by Hillebrand & Wolff [21] 
did not. Bruenn et al. [22]; also using the softest LS EOS variant and quasi-relativistic grav- 
ity, found strong explosions setting in within ^^250 ms after bounce in progenitors with ZAMS 
masses of (12, 15, 20, and 25) Mq. Suwa et al. [33], using Newtonian gravity and the soft LS 
EOS variant, found early, but weak explosions in a IS-M© progenitor star. Ott et al. [23] and 
Burrows et al. [25[ I26j , on the other hand, who performed purely Newtonian calculations using 
the stiffcr H. Shen EOS, did not find neutrino-driven explosions in progenitors of 11.2 — 25 Mq. 

Given that Nature has a way to robustly (without fine tuning) explode at least a significant 
fraction, but probably most stars with ZAMS masses of '-^lO — 20 Mq [271 [28], the large range of 
differing and sometimes disagreeing results of 2D simulations is dissatisfactory, if not disturbing. 

There are essentially three possible ways out: 
(1) The neutrino mechanism, while getting much closer to being viable in 2D than in ID, may 
still not be reaching its full efficacy. In 3D, an additional fluid motion degree of freedom is 
available and the nature of turbulence changefj^ This may allow accreting material to stay 
even longer in the region of net heating, resulting in a greater heating efficiency and, thus, 
potentially make the neutrino mechanism robust. Results to this effect have been obtained by 

^Provided that the turbulent cascade is resolved, turbulent power will cascade towards small scales in 3D 
while it cascades to large scales in 2D, which is unphysical. 
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Nordhaus et al. [29] who performed ID, 2D, and 3D calculations with parameterized neutrino 
heating and cooling using spherical Newtonian gravity and the H. Shen EOS. This work con- 
firmed the results of [13] for the ID— >2D case and found another big increase in efficacy when 
going from 2D to 3D. However, a similar parameterized study, carried out by Hanke et al. |16| , 
found no significant difference between 2D and 3D. The debate thus remains open and more 
work will be needed before the final word on the neutrino mechanism can be spoken. For 
this, fully self-consistent 3D simulations with reliable energy-dependent neutrino transport will 
be necessary. The first steps towards such self-consistent 3D models have already been taken 
[5U1 [5T] and their results, while not definite, are encouraging. 

(2) If dimensionality is not the key to robust neutrino-driven core-collapse supernova ex- 
plosions, then could there be physics missing from current ID and 2D simulations that, once 
included, could render 2D, or perhaps even ID, neutrino-driven explosions robust? A key ex- 
ample for this are self-induced (by p-p scattering) collective neutrino oscillations and we will 
discuss their potential effect on the neutrino mechanism in S|3] 

(3) If 3D and/or new physics cannot save the neutrino mechanisms, alternatives must be 
sought. Potential ones include the magnetorotational mechanism (e.g., [32]), the acoustic mech- 
anism [251 1261 155] . and the phase-transition-induced mechanisms f34|. The magnetorotational 
mechanism requires very rapid rotation in combination with non-linear magnetic field ampli- 
fication after bounce by the magnetorotational instability (e.g., [311 133 )■ Pulsar birth spin 
estimates [3S] and stellar evolution calculations that take into account magnetic fields (e.g., 
[37]) suggest that it may be active in no more than ~1% of massive stars that produce very 
energetic explosions and are related to the hyper-energetic core-collapse supernova explosions 
associated with a growing number of long gamma-ray bursts [5HHiO] . 

The acoustic mechanism^ proposed by |251 126] , relies on the excitation of protoneutron star 
pulsations by turbulence and SASI-modulated accretion downstreams. These pulsations reach 
large amplitudes at 600—1000 ms after bounce and damp via the emission of strong sound waves 
that steepen to secondary shocks as they propagate down the radial density gradient in the 
region behind the stalled shock. They dissipate and heat the postshock region, robustly leading 
to explosions, which, however, tend to be weak and occur late. This mechanism has not been 
confirmed by other groups, has been studied only in 2D simulations, and, most importantly, [41] 
have shown via non-linear perturbation theory that a parametric instability between the main 
mode of pulsation and abound higher-order modes, which are not resolved by the numerical 
models of 25, 26 , is likely to limit the mode amplitude to dynamically negligible magnitudes. 

The phase-transition-induced mechanism (e.g., |341I42] ) requires a hadron-quark phase tran- 
sition occurring within the first few 100 ms after core bounce (hence, at moderate protoneutron- 
star central densities). This phase transition leads to an intermittent softening of the EOS, a 
short collapse phase followed by a second bounce launching a secondary shock wave that runs 
into the stalled shock and launches an explosion even in spherical symmetry. However, the 
needed early onset of the phase transition requires fine tuning of the quark EOS and leads to 
maximum cold neutron star gravitational masses inconsistent with observations [31 143) . 

In the remainder of this contribution to the proceedings of the HAmburg Neutrinos from 
Supernova Explosions 2011 (HAi^SE 2011) conference, we discuss, in ^J2j new boundary condi- 
tions of core-collapse supernova theory set by neutron star mass and radius constraints, and, 
in fj3| we summarize the recent rapid progress made by studies considering the potential effect 
of collective neutrino oscillations on the core-collapse supernova mechanism. In ^|4] we cricially 
summarize our discussion and highlight the new frontiers of core-collapse supernova theory. 
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Figure 1: Mass-radius relations for 10 publically available finite-temperature EOS along with 
several constraints. The EOS are taken from [181 144H48] and the Tolman-Oppenheimer-Volkoff 
equation is solved with T = 0.1 MeV and neutrino-less /3-equilibrium imposed. The family of 
LS EOS is based on the compressible liquid-droplet model [TH] while all other EOS are based 
on relativistic mean field theory. The nuclear theory constraints of Hebeler et al. [JS] assume a 
maximum mass greater than 2 Mq and do not take into account a crust (which would increase 
the radius by ~400m). EOS that do not support a mass of at least 1.97 ± 0.04 are ruled 
out [21 133] . Ozel et al. [SD] analyzed three accreting and bursting neutron star systems and 
derived mass-radius regions shown in green. Steiner et al. |51j performed a combined anaylsis 
of six accreting neutron star systems, shown are 1-cr and 2-a results in blue. 



2 New Constraints on the Supernova Equation of State 

An important ingredient in any core-collapse supernova model is the nuclear EOS. It provides 
the crucial closure for the set of (magneto-)hydrodynamics equations used to describe the 
evolution of the collapsing stellar fluid and strongly influences the structure of the protoneutron 
star and the thermodynamics of the overall problem. Nuclear statistical equilibrium (NSE) 
prevails above temperatures of ~0.5 MeV, which corresponds to densities above ~10^ gcm~'^ in 
the core-collapse supernova problem. In this regime, the EOS is derived from the Helmholtz free 
energy and thus is expressed as a function of density p, temperature T, and electron fraction 
Ye- The NSE part of the core-collapse supernova EOS must cover tremendous ranges of density 
(10^ — 10^^ gcm^'^), temperature (0.5 — 100 MeV), and electron fraction (0—^0.6). Constraints 
from experimental nuclear physics on the nuclear EOS are few and generally limited to only 
small regions of the needed (p, T, Ye) space (see the discussions in [47ll52]). 

A stringent constraint on the nuclear EOS is set by precision mass measurements of neutron 
stars in binary systems. The 2-M0 ([1.97 ±0.04] Mq) neutron star of Demorest et al. [H] rules 
out a large range of soft hadronic, mixed hadronic-exotic, and strange-quark matter EOS [3ll53j. 
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Recently, Hebeler et al. [49] have carried out chiral effective field theory calculations of 
neutron-rich matter below nuclear saturation density, strongly constraining the P{p) relation- 
ship in this regime. They derived a radius constraint for a 1A-Mq neutron star of 10.5 km < 
R < 13.3 km (these numbers would be shifted up by ^400 m if a detailed crust treatment was 
included) by requiring that all EOS support neutron stars with mass > 2 Mq and pass through 
the P{p) range allowed by their calculations. 

Steiner et al. |51| and Ozel et al. 150] analyzed observations from accreting and bursting 
neutron stars to obtain neutron star mass-radius constraints. Such observations and their 
interpretations should be taken with a grain of salt, since large systematic uncertainties are 
attached to the models that are required to infer mass and radius and to the assumptions made 
in their statistical analysis. For example, [51) and [50j . starting with different assumptions, 
derive rather different 2-cr mass-radius constraints from the same set of sources. 

In Fig. [l] we contrast the various observational constraints on the neutron star mass and 
radius with a range of EOS used in core-collapse supernova modeling. The LS family of EOS is 
based on the compressible liquid droplet model [THj, while all other EOS (drawn from [^1^148) ) 
are based on relativistic mean field (RMF) theory. The details of the M — R curves depend 
on multiple EOS parameters such as nuclear incompressibility, symmetry energy and their 
derivatives and we must refer the reader to [47] and to the primary EOS references for details 
for each EOS. Fig. [ij shows that none of the current set of available EOS allow for a 2-Mq 
neutron star while at the same time being consistent with the current mass-radius constraints 
from observations. The crux is that the EOS needs to be sufficiently stiff to support 2-Mq 
neutron stars and at the same time sufficiently soft to make neutron stars with moderate radii 
in the canonical mass range. This balance appears to be difficult to realize. The stiff set of 
RMF EOS produce systematically too large neutron stars. The soft compressible liquid-droplet 
LS180 EOS [T51 agrees well with the mass-radius constraints, but is ruled out by its failure to 
support a 2-Mq neutron star. Closest to satisfying all constraints are the LS220 EOS of [TH] 
and the yet unpublished HSDD2 EOS of [IS] based on the RMF model of [M]- 

The stiffness of the nuclear EOS at high and intermediate densities has important conse- 
quences for the postbounce evolution of core-collapse supernovae. In simple terms: the stiffer 
the EOS, the more extended the protoneutron star and the larger the radius and the lower the 
matter temperature at which neutrinos decouple from the protoneutron star matter. Assuming 
a Fermi-Dirac spectrum with zero degeneracy, the mean-squared energy of the emitted neutri- 
nos is approximately given by (e^) ~ 21T^ [55] . where is the matter temperature (in units of 
MeV) at the neutrinosphere (where the optical depth r « 2/3). Hence, a softer EOS will lead 
to systematically harder neutrino spectra than a stiffer EOS (as born out by the simulations of 
[17]V Since the charged-current neutrino heating rate Q+ scales oc (e^), a soft EOS leads to a 
higher neutrino heating efficiency than a stiff EOS. This is at least part of the explanation why 
some published 2D simulations using the soft, now ruled-out LS180 EOS have shown neutrino- 
driven explosions [TT] [HI [53] while simulations with stiffer EOS have generally failed to yield 
such explosions in stars more massive than '-^ll Mq [T71 [TOl [25] . 

3 New Physics: Collective Neutrino Oscillations 

Neutrinos and antineutrinos of all three flavors are produced in core-collapse supernovae and can 
oscillate from one flavor to another. Ve and are made and interact via charged-current and 
neutral-current interactions, while and Vr and their antineutrinos experience only neutral- 
current processes, since no muons or tauons are present in the core-collapse supernova environ- 
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ment. Hence, their interaction cross sections are very similar and one generally lumps them 
together as = {i/^, v^) and = {f^^, Vr}. 

The oscillations between and or v^, and are driven by their mass differences, forward 
scattering off background electrons, and forward scattering off other neutrinos and antineutri- 
nos. These limiting regimes are called neutrino oscillations in vacuum |56| . matter-enhanced 
oscillations through the MSW effect [57], and collective oscillations [58], respectively. Quanti- 
tatively, the nature of neutrino flavor conversions depends on an interplay of vacuum neutrino 
oscillation frequency = /S.m? / {2E) with the matter potential A = \f2Gpne due to back- 
ground electrons (where is the electron number density) and with the collective neutrino 
potential ^ ~ y/^Gp'y^ — cos0)nj^+p generated by the neutrinos themselves (where n^x+p is the 
neutrino and antineutrino number density). In a typical core-collapse supernova environment, 
the matter potential falls off with radius as Uf. oc whereas the collective potential falls off 

faster with n^^j^{l — cos 6*) cx l/r"*. So, when the neutrinos travel outward from the core, they 
generally first experience collective effects and then matter effects, which may be modified by 
shock wave effects |59j. After they leave the star, the mass eigenstates travel independently 
and are detected on Earth as an incoherent superposition. There can be distinctive effects due 
to additional conversions during propagation inside the Earth (e.g., |60j). 

3.1 Collective Oscillations due to u-u Interactions 

The neutrino density creates a potential that is not flavor diagonal '58J; n^^ rip are density matri- 
ces in flavor space and depend on the flavor composition of the entire neutrino ensemble! Flavor 
evolution of such dense neutrino gases |61j can be understood to good accuracy without consid- 
ering many-particle effects p51 . Calculations in spherical symmetry showed that the collective 
oscillations can affect neutrino flavor conversions substantially [S31 |M| . The main features ob- 
served were large flavor conversions for inverted hierarchy (neutrinos masses mi,m2 > 7713), 
and a surprisingly weak dependence on the mixing angle and the matter density. 

These features can be understood analytically. A dense gas of neutrinos displays collective 
flavor conversion [HS], i.e., the flavor oscillations of all neutrinos and antineutrinos become cou- 
pled to each other and all of them undergo flavor conversion together. Neutrinos of all energies 
oscillate almost in phase, through synchronized |66j / parametricallv resonant [67] /bipolar os- 
cillations [5H1 [M] • The effect of the bipolar oscillations with a decreasing collective potential 
/X is a partial or complete swap of the energy spectra of two neutrino flavors |701 171j . The 
"1 — cos 9" structure of weak interactions can give rise to a dependence of flavor evolution on 
the neutrino emission angle [Mj or even flavor decoherence, i.e., neutrinos acquire uncorrelated 
phases, and the neutrino fluxes for all flavors become almost identical [72]. For a realistic ex- 
cess of i^e, compared to fluxes, such angle-dependent effects are likely to be small [751 174j . 
Even non-spherical source geometries can often be captured by an effective single-angle approx- 
imation [75\ in the coherent regime. While most of these results were obtained for neutrino 
oscillations between two flavors, it was shown that with three flavors one can usually treat the 
oscillation problem by factorizing it into simpler two-flavor oscillation problems, since the mass- 
squared differences between the mass eigenstates obey Am^2 ^ l^'Tiisj and the mixing angle 
^13 ^ 1 [ZSjj and the previous results are easily generalized. Effects of potential CP violation 
are expected to be small with realistic differences between /i and r neutrino fluxes |77| . On 
the other hand, similar realistic departures are sufficient to trigger collective effects even for a 
vanishing mixing angle [78l [79] . 
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3.2 Results obtained with Core-Collapse Supernova Toy-Models 



Although the inherent nonhnearity and the presence of multi-angle effects make the analysis 
rather complicated, the final outcome for the neutrino fluxes turns out to be rather straightfor- 
ward, at least in the spherically symmetric scenario. Synchronized oscillations with a frequency 
(w) take place just outside the neutrinosphere at r ~ 10 — 40 km. These cause no significant 
flavor conversions since the mixing angle, which determines the extent of flavor conversion, is 
highly suppressed by the large matter potential due to the high electron density in these in- 
ner regions |801 181j . A known exception occurs for the burst phase in low-mass progenitor 
stars that have a very steep density profile [55]. In such a situation, neutrinos of all energies 
undergo MSW resonances before collective effects become negligible [531 [51] . At larger radii, 
r ~ 40— 100 km, bipolar or pendular oscillations i>e o with a higher frequency y^2ujfi follow. 
These oscillations are instability-driven and thus depend logarithmically |68j on the mixing an- 
gle, occurring where the fluxes for the two flavors are very similar [71]. As /x decreases, so that 
(a;) ^ fi, neutrinos near this instability may relax to the lower neutrino mass (energy) state. As 
a result, one finds one or more spectral swaps demarcated by sharp discontinuities or "spectral 
splits" in the oscillated flux. 

These simple explanations do not take into account the fact that neutrinos are emitted 
at different angles from the neutrinosphere. As a result, radial neutrinos take a shorter path 
(while tangentially emitted neutrinos take a longer path), and thus experience less (more) 
background potentials from the electrons and from other neutrinos leading to an emission 
angle-dependent flavor evolution. These sort of effects are called multi-angle effects, and can 
suppress or delay flavor conversions either through multi-angle matter effects [85] . or through 
multi-angle neutrino-neutrino interactions themselves |86j . 



3.3 Results obtained with more realistic Models 



The observations outlined in the previous sect ion [3 . 2 1 were mostly based on toy models of core- 
collapse supernova neutrino fluxes and background densities. Recently, several groups have 
tried to perform semi-realistic calculations of the oscillation physics, by injecting the output 
neutrino fluxes from supernova simulations into oscillation calculations |87H89| . Interesting 
results have also been obtained by performing a linear stability analysis of the equations used 
for calculating the flavor conversion, with the initial conditions taken from simulations j90[ I91j . 

The simple picture given in |3.2| has therefore undergone further changes. Firstly, it has 
been recognized that matter effects suppress collective oscillations even in the bipolar regime 



through multi-angle effects as explained before [571 [55] (see also ^3.4). This is most effective 
in the pre-explosion accretion phase, when the matter density is large in the region behind the 
stalled shock. In this case, it appears that one can simply ignore collective effects and only 
include the MSW effects which take place at larger radii. Of course, the result depends on 
details of the matter density and ratios of neutrino fluxes. In particular, for fluxes that are 
either highly symmetric in neutrinos and antineutrinos j89| , or include flavor dependent angular 
emission that leads to an angular instability [HH |92] , one flnds the matter suppression to be less 
effective. Secondly, in the cooling phase of the explosion, v^jv^ fluxes may be larger than Ve/^e 
fluxes. This can lead to additional instabilities which cause multiple spectral splits [7T]. These 
features survive multi-angle effects in general, and with the inclusion of three-flavor effects can 
lead to a rich and complex phenomenology [SSI [SI] ■ 

The understanding of collective neutrino oscillations is still evolving, and we expect that 



Proceedings ofHAvSE 2011 



7 



1000 



100- 



o 



10; 



1 - 



0.1 



0.01 



0.001 



S11.2WHW02 - expected (single-angle) 
si 1.2WHW02 - expected (multi-angle) 
S11.2WHW02 - < ig^„ (hypothetical) 
S15WH07 - expected (single-angle) 
S15WH07 - expected (multi-angle) 
S15WH07 - r,„j < r (hypothetical) 




0.15 0.2 0.25 

t ~ tbounce [s] 



0.4 



Figure 2: Time evolution (as a function of time after core bounce) of the potential percentage 
increase in the heating rate due to collective neutrino oscillations based on our recent simulations 
[HI], in which we considered 11.2-M0 and 15-Mq progenitors. The dashed lines assume the naive 
case of complete conversion already below the gain radius where heating begins to dominate 
over cooling. This is the case assumed by by Suwa et al. [HS]) and leads to an enhancement 
of up to 100%. In the our detailed oscillation calculations, conversion does not occur before 
the gain radius and our more realistic estimate of the heating enhancement is much lower and 
shown in solid lines. The points, blue squares for the 11.2-M0 model and red circles for the 15- 
Mq model, represent our estimate of the heating enhancement if multi-angle oscillation effects 
are included, which further increase the radii at which collective oscillations occur and thus 
decrease the heating enhancement even further, in general agreement with |87l |88j . This figure 
corresponds to Fig. 7 of [89] . 



more accurate numerical calculations and improved analytical understanding will yield new 
surprises and insights into the existing results we have summarized here. 

3.4 Effect of Collective Oscillations on Neutrino-driven Explosions 

From the core-collapse supernova theory point of view, the most intriguing result of collective 
oscillations is the almost complete exchange of Ve and Vx and and spectra in the inverted 
mass hierarchy. The v,j. and are emitted by thermal processes deep inside the core and 
their spectra are much harder than those of their electron-flavor counterparts. Due to the e^- 
dependence of the charged-current absorption cross section, a swap of Vx/^x and Ve/i^e spectra 
could dramatically enhance neutrino heating and may be the crucial ingredient missing in core- 
collapse supernova models, provided that the oscillations occur at sufflciently small radii to 
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have an effect in the region behind the shock. To our knowledge, this point, in the context of 
collective oscillations, was first made by one of us 

Suwa et al. jHl] recently performed a set of ID and 2D core-collapse supernova simulations 
in which they considered ad-hoc spectral swaps above 9 MeV for neutrinos and antineutrinos 
occurring at a fixed radius of 100 km, which is close to the gain radius (where heating begins 
to dominate over cooling) in their simulations. They considered a range of progenitor models 
and found that the heating enhancement by collective oscillations can indeed turn duds into 
explosions. This result was corroborated in a semi-analytic study by Pejcha et al. ^97J in which 
the authors also considered different radii for the oscillations to become effective. 

Chakraborty et al. [87|[88] carried out the first multi-angle single-energy neutrino oscillations 
based on realistic neutrino radiation fields from ID core-collapse supernova simulations. They 
discovered that the rather high matter density between protoneutron star and stalled shock 
strongly suppresses collective neutrino oscillations in the pre-explosion phase when multi-angle 
effects are taken into account. Hence, the authors excluded any impact of collective oscillations 
on neutrino heating. 

In Dasgupta et al. [89], we carried out single-angle multi-energy and multi-angle single- 
energy oscillation calculations based on neutrino radiation fields from 2D core-collapse super- 
nova simulations performed with the VULCAN/2D code [25]. In 2D, convection and SASI lead 
to complicated flow patterns and large-scale shock excursions not present in ID simulations. 
Even in our single-angle calculations and in the most optimistic case, we find that collective 
oscillations do not set in at radii sufficiently deep in the heating region to have a significant 
effect on neutrino heating. When including multi-angle effects, we also observe a suppression 
of collective oscillations, though not at the level argued for by [57] [55] i who made different 
assumptions about the angular distribution of the neutrino radiation fields emitted from the 
neutrinosphere (ours are based on the angle-dependent neutrino transport results of [2^). 

As depicted by Fig. [2] we find that the heating enhancement due to collective oscillations, if 
present at all, stays below ~0.1% at all times in both considered progenitor models when oscil- 
lation radii from full oscillation calculations are taken into account. This shows, in agreement 
with |871 188j , that the strong positive effect on the neutrino mechanism reported by Suwa et 
al. [nn] is artificial and due primarily to their ad-hoc choice of a small oscillation radius. 



3.5 Collective Oscillations after the Onset of Explosion 

In Dasgupta et al. [53], we studied the suppression of collective oscillations by multi-angle 
effects at high matter density using multi-angle single-energy oscillation calculations for select 
simulation snapshots of the pre-explosion phase in 11.2-Af0 and Ih-M^ progenitors. We also, 
in a more heuristic approach, studied the potential for suppression of collective oscillations 
by comparing the MSW potential \{r) with the expression Ama = 2-\/2G_F<I'i/,p(i?^^/r^)J^_, 
where ^ is the neutrino number density at the Vf. neutrino sphere radius R^^ and J-"_ = 
{^u^ — *i'pe)/(*J'i'c + '^Pe + 4<I>,y^) Is thc rclativc lepton asymmetry of the neutrinos. If A ^ Ama, 
collective oscillations are suppressed j85] . 

In [89], we compared A and Ama at various pre-explosion times, radii, and spatial angular 
directions in our simulations using 11.2-M0 and 15-M0 progenitors. Based on this, we con- 
cluded, in agreement with [571 [55], that suppression of collective oscillations is likely highly 
relevant in the pre-explosion phase and must be carefully studied even in relatively low-mass 
progenitors with steep density profiles such as the 11.2-Af0 progenitor model. 

The situation after the onset of explosion, however, may be quite different: The explosion 
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Figure 3: The MSW potential A(r) along various directions (solid lines, 10 rays equally spaced 
in cos[0iat]), in comparison to the minimum A(r) needed for multi-angle suppression, Ama = 
2y/2G F^u,v{R^^/ r^)^-- Dot-dashed-dashed, dot-dot-dashed, and dashed lines, indicated Ama 
taken along the North pole (NP), equator (Eq), and South pole (SP), respectively. The steep 
rise in the A(r) profiles at 250 ms and 350 ms occurring around r ^ 1000 km and r ~ 2000 km, 
respectively, is the location of the shock. At 450 ms the shock is close to 3000 km. 



rarefies the region behind the expanding shock and shuts off the large I'e/i'e accretion lumi- 
nosity, changing the neutrino flux asymmetry. Extending our previous results to the explosion 
phase, we have repeated our simulations for the 11.2-Mq progenitor, but included an additional 
neutrino heating term in order to drive an early explosion. The heating term is equivalent to 
the prescription used in [T3] with L^^ = = 0.5 x 10^^ ergs/s. 

In Figure |3j we compare the MSW potential A(r) along multiple angular directions with 
Ama- As the explosion clears out the region behind the shock, the MSW potential decreases in 
strength. In this model, within a few 100 ms of the onset of the explosion, the MSW potential 
becomes comparable to AmAj which indicates that the suppression is lifted and collective oscil- 
lations may now occur at radii as small as ~200km. At this point, the core-collapse supernova 
explosion has already been launched, but collective neutrino oscillations may still affect the 
evolution and various observable features, for example, via the neutrino-driven wind from the 
protoneutron star and r-process nucleosynthesis (98, ,99j . 

4 Summary and Outlook 

In this contribution to the proceedings of the HAmburg Neutrinos from Supernova Explosions 
(HAi/SE) 2011 conference, we have summarized the recent rapid progress in various aspects of 
core-collapse supernova theory. While 2D simulations continue to be perfected IMl [TIM 1101 j . 
self-consistent 3D simulations with energy-dependent neutrino radiation hydrodynamics are 
now the frontier of core-collapse supernova modeling [30] and are made possible by the first 
generation of petascale supercomputers. General relativity is also beginning to be included 
in 2D (1UD[4102] and 3D simulations [1031 1104j , which will eventually allow for first-principles 
studies of multi-D black hole formation and the relationship between massive star collapse 
and long gamma-ray bursts. Also, open-source codes and microphysics inputs (EOS, neutrino 
opacities) are gaining traction [48l ll05H112j . They allow for code verification and physics 
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benchmarking and are lowering the technological hurdle for new groups with new ideas trying 
to enter core-collapse supernova modeling. 

After ten years with little activity, improved modeling capabilities, faster computers, and 
the discovery of the 2-Mq neutron star, have spawned a flurry of activity in the nuclear EOS 
community, which has already resulted in multiple new finite-tempcraturc EOS for core-collapse 
supernova modeling [iiHlT] . 

The realization that collective neutrino oscillations may occur in the core-collapse supernova 
environment [MIIMIIM] has led to a plethora of work since ^2005. As we have outlined in this 
article, the current state of affairs is that collective oscilations are unlikely to be dynamically 
relevant in driving the explosion, but their effects are crucial in predicting and understanding 
the neutrino signal that will be seen in detectors from the next nearby core collapse event. 
The current frontier of oscillation calculations in the core-collapse supernova context is marked 
by detailed multi-energy multi-angle calculations that take their input spectra and angular 
distributions from core-collapse supernova models. Significant progress towards this has re- 
cently been made [87^f89l I113j , but more will be needed to assess the potentially strong impact 
of neutrino-matter interactions and only partially-decoupled neutrino radiation fields in the 
oscillation regime. 

The broad range of current and near-future advances in theory will be matched and tested 
by observations of the next galactic (or Magellanic-cloud) core-collapse supernova. This event 
will most likely be observed in electromagnetic waves, neutrinos, and, with the upcoming ad- 
vanced generation of gravitational-wave observatories, for the first time also in gravitational 
waves. Gravitational waves carry dynamical information on the intricate multi-D processes oc- 
curing in the supernova core [1141 1115] and will complement the structural and thermodynamic 
information carried by neutrinos. Together, neutrinos and gravitational waves may finally shed 
observational light on the details of the core-collapse supernova mechanism. 
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